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In this paper we perform stability analysis for exponential solutions in Einstein-Gauss- 
Bonnet and cubic Lovelock gravity. We report our findings, provide areas on parameters 
space and discuss familiarities and differences between cases. Analysis suggests that only 
several cases out of numerous found solutions could be called stable. In particular, cases with 
three-dimensional isotropic subspace which could give rise to successful compactification are 
diminished to one general case and one additional partial solution in the cubic Lovelock case. 
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I. INTRODUCTION 

Recent developments in particle physics, especially experiments on Large Hadron Collider, bring one 
closer and closer to the lowest energy scale when one can start to probe extra dimensions. Despite the fact 
that by now only Kaluza-Klein theories are somehow tested, it is just a matter of time before we will be able 
to test many more extradimensional theories with particle physics. In that regard alternative tests are more 
then welcommed and cosmology could be of use here. One of the simplest implementations of cosmology in 
this question is search for cosmological models which allow compactification. 

One of the first attempts to find an exact static solution with metric being a cross product of a (3-1-1)- 
dimensional manifold and a constant curvature “inner space”, also known as “spontaneous compactification”, 
were done in [1], but with (3-|-l)-dimensional manifold being actually Minkowski (the generalization for a 
constant curvature Lorentzian manifold was done in [2]). In the cosmological context it could be useful 
to consider Priedman-Robertson-Walker metric for (3-|-l)-dimensional section; this situation with constant¬ 
sized extra dimensions was considered in [3] . There it was explicitly demonstrated that to have more realistic 
model one needs to consider the dynamical evolution of the extra dimensional scale factor as well. In the 
context of exact solutions such an attempt was done in [4] where both (3-1-1)- and extra dimensional scale 
factors were exponential functions. Solutions with exponentially increasing (3-|-l)-dimensional scale factor 
and exponentially shrinking extra dimensional scale factor were described. 

More recent analysis focuses on properties of black holes in Gauss-Bonnet [5, 6] and Lovelock [7, 8] 
gravities, features of gravitational collapse in these theories [9-11], general features of spherical-symmetric 
solutions [12] and many others. Of recent attempts to build a successful compactification particularly rel¬ 
evant are [13] where the dynamical compactification of (5-1-1) Einstein-Gauss-Bonnet (EGB) model was 
considered, [14, 15], with different metric ansatz for scale factors corresponding to (3-1-1)- and extra dimen- 


sional parts, and [16, 17] where general (e.g. without any ansatz) scale factors and curved manifolds were 
considered. 

In [2] the structure of the equations of motion for Lovelock theories for various types of solutions has 
been studied. It was stressed that the Lambda term in the action is actually not a cosmological constant as 
it does not give the curvature scale of a maximally symmetric manifold. In the same paper the equations 
of motion for compactification with both time dependent scale factors were written for arbitrary Lovelock 
order in the special case that both factors are flat. The results of [2] were reanalyzed for the special case of 
10 space-time dimensions in [18]. In [19] the existence of dynamical compactification solutions was studied 
with the use of Hamiltonian formalism. 

Usually when dealing with cosmological solutions in EGB or more general Lovelock gravity [20] one 
imposes a certain ansatz on the metric. Two most used (and so well-studied) are power-law and exponential 
ansatz. The former of them could be linked to Friedman (or Kasner) stage while the latter - to inflation. 
Power-law solutions were intensively studied some time ago [2, 21] and recently [22-26] which leads to almost 
complete their description (see also [27] for useful comments regarding physical branches of the solutions). 
Exponential solutions, on the other hand, for some reason are less studied but due to their “exponentiality” 
could compactify extra dimensions much faster and more reliably. Our first study of exponential solutions [28] 
demonstrate their potential and so we studied exponential solutions in EGB gravity full-scale. We described 
models with both variable [29] and constant [30] volume and developed general solution-building scheme 
for EGB; recently [31] this scheme was generalized for general Lovelock gravity of any order and in any 
dimensions as well. 

In some sense this paper is a logical continuation of [29, 31] - in there we found analytically exponential 
solutions in Einstein-Gauss-Bonnet and cubic Lovelock gravity and now we want to perform stability analysis 
of the solutions found. This is important since just finding analytical solutions often is not enough - if the 
solution found is unstable, it cannot describe physical phenomena it is intended to describe. In our case it is 
compactihcation - indeed, despite the fact that in [29, 31] we described all possible spatial splittings, special 
attention was paid to cases with three-dimensional isotropic subspace, which, being expanding while all 
remaining dimensions are contracting, give us good example of compactification scheme at work. In [29, 31] 
we found a number of solutions with three-dimensional isotropic subspace - now it is time to check their 
stability to see how much of them are stable. 

For stability analysis we use four cases - (4-|-l), (5-1-1), (6-1-1) and (7-1-1) dimensions. The reason for 
this choice is quite simple - as we stated, we analyze the stability of the solutions found in [29, 31] and 
they are found in these four different dimensions. And the reason for considering these four in first place 
is also simple - (4-|-l) and (5-1-1) are the lowest two dimensions where Gauss-Bonnet gravity plays a role 
while (6-|-l) and (7-1-1) are two lowest dimensions where cubic Lovelock plays a role. So that we not only 
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consider solutions in two lowest dimensions, but also compare them in (6+1) and (7+1) for EGB and cubic 
Lovelock cases to see the difference cubic Lovelock term brings. 

The structure of the manuscript is as follows: first we consider Gauss-Bonnet case and study stability 
for (4+1)-, (5+1), (6+1)- and (7+l)-dimensional cases consequently, followed by a brief conclusion on the 
considered cases. Then we deal with solutions with cubic Lovelock term in (6+1) and (7+1) dimensions, 
also followed by a brief conclusion. Finally, we discuss our results, point out familiarities and differences 
between cases and build directions for further research. 

II. GAUSS-BONNET SOLUTIONS 

In this section we study stability of the Gauss-Bonnet exponential solutions found in [29, 31]. We 
will separately consider cases with different number of spatial dimensions and briefly comment each case. 
As [29, 31] were devoted to study of vacuum and A-term solutions only, the equations of motion are written 
in a way to fit only these two cases - with A-term as a source; in that case vacuum solutions are obtained 
with A = 0. Full Einstein-Gauss-Bonnet system (i.e. without exponential ansatz) reads: jth dynamical and 
constraint equations 


2ij:{H, + Hf)+ E HiHk] +8a[j:iHi + Hf) E HkHi + 3 E HkHiHmHn] - A = 0 

2 E HiHj + 24q; E HkHiHmHn = A, 

i>j k>l>m>n 

( 1 ) 

where Hi = a is Gauss-Bonnet coupling and we put Einstein-Hilbert coupling to unity. 

We perturb full system (1) around exact exponential solutions to find the regions of variables and pa¬ 
rameters where they are stable. As we perturb solution as Hi ^ Hi + 5Hi and Hi 5Hi (since exact 
exponential solution has H = 0), and keeping in mind that equations of motion are first order, the criterium 
for stability would be 


SHi 

6Hi 


< 0 . 


So we derive 5Hi/6Hi and find regions on variables and parameters space which satisfy (2). 


( 2 ) 


A. 4+1 

According to [29], (4+l)-dimensional case has solutions with three different spatial splittings: isotropic 
(4 + 0) = {H, H,H,H}, (3 + 1) = {H, H, H, h} and (2 + 2) = {H, H, h, h}. First of them, isotropic, is the 
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Figure 1. Stability regions on the {H, a) plane for Gauss-Bonnet solutions. In (a) panel we presented typical behavior 
for {{D — 2) + 2) spatial splitting. The situation in (a) panel exactly correspond to (4-1-1 )-dimensional (2-1-2) = 
{H, H, h, h} case (see Eq. (4)), but other cases with {{D — 2) -|- 2) spatial splitting - (3 -|- 2) = {H, H, H, h, h} case 
from (5-1-1) (see (6)), (4 -|- 2) = {H, H, H, H, h, h} case from (6-1-1) (see (7)) and (5 -I- 2) = {H, H, H, iJ, H, h, h} case 
from (7-1-1) (see (9)) - have similar behavior - they meet this picture up to a factor. In (b) panel we have typical 
behavior for ((D — 3) -I- 3) spatial splitting. Again, exactly this figure corresponds to (3 -I- 3) = {H, H, H, h, h, h} case 
from (6-1-1) dimensions (see Eq. (8)), but up to a factor the situation with (4 -|- 3) = {H, H, H, H, h, h, h} case from 
(7-1-1) dimensions (see Eq. (10)) could also be described by this figure (see text for details). 


simplest one - with only one variable we easily get 6Hi = —AHSHi (and the perturbed equations are subject 
to cyclic indices permutations) so that the solution is stable as long as > 0; this is true for both vacuum 
and A-term solutions. 

Second to consider, (2 -|- 2) = {H, H, h, h} is “symmetric” in a way that interchanging H and h does not 
change equations of motion. Solutions of the perturbed equations read 


= (3) 

so that stability is reached as long as {H -|- h) > 0. Remembering that h = —l/(4af7), we can find stability 
regions on the {H, a) space - they read 


4ai/2 _ I 

AaH 


> 0 


( 4 ) 


and are presented in Fig. 1(a). 

Final solution, (3 -|- 1) = {H,H,H,h}, gives rise to an interesting situation: solution reads [29] /i S R, 
and so when we perturb it the perturbation which corresponds to h appears to be undefined - 5H4^ could be 
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eliminated from the resulting system of the perturbed equations, making the perturbation in that direction 
undefined. Formally since perturbation is undefined, we cannot call it “unstable”, but on the other hand it 
is definitely not stable either, which in a sense makes it in some sort of neutral balance. In a sense this is 
expected - indeed, since h is unconstrained, one has all grounds to question physical validity of the solution, 
so it is just expected that the solution is not stable. 

To conclude, two out of three (4+l)-dimensional solutions have stability regions on the {H, a) plane: 
isotropic solution is stable iff > 0 while (2 + 2) solution is stable iff {H + h) >0; considering the solution 
for h one can find stability regions on the {H, a) space - they are presented in Fig. 1(a). Final solution to 
consider, (3 + 1), the only one which could give a rise to viable compactification scheme, appears to have 
unconstrained perturbation (due to lack of constraint on /i E R), so that we cannot put it into “stable” class. 
So that further we call these solutions “unstable” but only in a sense that we cannot definitely call them 
“stable”. 


B. 5+1 

(5+l)-dimensional case has, according to [29], solutions with four different spatial splitting: isotropic 
(5 + 0) = (4 + 1) = {H,H,H,H,h} case, (2-2 + 1) = {H, H,-H,-H, h} case and 

(3 + 2) = {H, H, H, h, h} case. 

First of them - isotropic one (5 + 0) = {H, H, H, H, H} - is quite similar to the (4+l)-dimensional case. 
Again, perturbed equation reduced to dHi = —5H6Hi and so similar to (4+l)-dimensional case the solutions 
(both vacuum and A-term) are stable iff > 0. 

The second case - (4 + 1) = {H, H, H, H, h} - is also similar to the analogue from (4+1) dimensions - 
due to unconstrained /i E R we have excluded from the system of perturbed equations and so unstable 
with respect to the perturbations in the direction. 

The next case - (2 — 2 + 1) = {H, H, —H, —H, h} - is similar to (4 + 1) = {H, H, H, H, h} in a sense that 
we again have /i E R and so the solution is unstable in direction. 

The last remaining splitting (3 + 2) = {H, H, H, h, h} has different solutions for A-term and vacuum 
cases. Despite of that the latter could be considered as a case of the former, and it has the solution for 
perturbations: 


6H, 


Qe 


4 aH J 


(5) 


Again, remembering that h = + l)/{8aH), we can draw regions of stability on the {H, a) space 

- they read 


5 


( 6 ) 


8aH^ — 1 
AaH 


> 0 ; 


one can see that up to a factor it behaves similar to (4) and so qualitatively it is the same as presented in 
Fig. 1(a) 

Final case, vacuum (3 + 2) = {H, H, H,h,h}, is a particular case from A-term case, it is obtained just 
with A = 0. In that case we have 192^^ — 122^^ + 4^ — 1 = 0 for ^ = aH^ and so ^ ^ 0.56276. Further, 
h = —{AaH'^ + l)/{8aH) ~ —0.72277 and the exponent in (5) become ~ —1.55677, making this solution 
stable for 77 > 0; on a plot it situated in the first quadrant inside stable area. 

This finalize our study of (5+l)-dimensional solutions stability. First two solutions are similar to the 
previous case - isotropic solution is stable as long as 77 > 0 while (4 + 1) = {77,77,77,77, h} is unstable in the 
/i-direction since h is unconstrained. We also have unphysical solution for (2 — 2 + 1) = {77,77, —77, —77, h} 
splitting case where h is also appears to be unconstrained in the general case. Finally, we found stability 
regions for (3 + 2) = {77, 77, 77, h, h} A-term case. Vacuum (3 + 2) = {77, 77, 77, h, K] solution appears to be 
stable iff 77 > 0; it is a particular case with A = 0. 


C. 6+1 

According to [31], (6+l)-dimensional case admits solutions in five spatial splittings: isotropic (6 + 0) = 
{77,77,77,77,77,77}, (5 + 1) = {77, 77, 77, 77, 77,/i} case, (3 + 2 + 1) = {H, H, H,h,h, z} case, (3 + 3) = 
{77,77,77, h, h, h} and (4 + 2) = {77,77,77,77, h, h} case. First two of them are similar to the previous cases 
- isotropic solution is stable iff 77 > 0 while (5 + 1) = {77,77,77,77,77, h} is unstable since /i € K.. Let us 
consider three remaining cases. 

First of them, (3 + 2 + 1) = {77,77,77, h, h, z} case, is unstable as 2 ; € R and through this it is unstable in 
2 ; direction. Last of them, (4 + 2) = {77,77,77,77, h, h} case, has the following stability regions for perturbed 
equations: 


12a772 - 1 
6aH 


> 0 ; 


(7) 


again, qualitatively it behaves similar to the case presented in Fig. 1(a). This splitting also admits vacuum 
solution which is found from general A-case as A —)• 0 limit. With ^ = h/H definition we have + 3^^ + 
6^ + 5 = 0 with solution ^ ~ —1.322. Further we have 77^ = — 1/(12q:(^ + 1)) which after reversion and 
substitution reads a ~ 1/(3.86477^). From (7) one can see that separation curves now are a = 1/(1277^) so 
that again vacuum curve intersect stable area in first quadrant and so vacuum solution is stable iff 77 >0. 
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Final case - (3 + 3) = {H, H, H, h, h, /i} - is a bit more interesting - its stability condition reads 


SHi 


Cie 




12aif2 _ 1 


t 


(8) 


and the corresponding areas are plotted in Fig. 1(b). The graph appears to be asymmetric - for H < 0 we 
have a > l/(12Ff^) while for Ff > 0 it is a > 1/(8FF^). One cannot miss that in (8) we have two signs while 
the graph in one - since we are interested in existence of stable solution, we put both branches together, so 
that at each point inside area at least one of branches is stable. 

This splitting also admits vacuum solution. Following the procedure from (4 + 2) = {H, H, H, H, h, h} 
case, we have — 3^ + 1 = 0 and its solutions are ^ = 3/2+ while solution for could be reverted to 

a = —l/(4iF^((^^+4^+l)). Substitution of ^ values found to expression for a gives us 5Hi/6Hi = 3/2(1+ \/5) 
so that branch is stable iff FF > 0 while is stable iff FF < 0. 

This finalize our study of (6+l)-dimensional solutions stability. As expected, isotropic case has the same 
behavior as in all previous cases while solutions with unconstrained exponents ((5 + 1) = {FF, FF, FF, FF, FF, K} 
and (3 + 2 + 1) = {FF, FF, FF, h, h, z} cases) have instabilities in the corresponding directions. First of final 
two cases - (4 + 2) = {FF, FF, FF, FF,/i, F} - has “typical” stability areas (see Fig. 1(a)) while the second 
one - (3 + 3) = {H, H, H,h,h,h} case - has its own (Fig. 1(b)). Vaccum counterpart of the (4 + 2) = 
{FF, FF, FF, FF, h, h} solution is stable iff FF > 0 while vacuum counterpart of (3 + 3) = {FF, FF, FF, h, h, h} case 
has two branches - one of them is also stable iff FF > 0 while the other one is stable only iff FF < 0. 


D. 7+1 

This case allows solutions with seven different spatial splittings: isotropic (7+0) = {FF, FF, FF, FF, FF, FF, FF}, 
(6 + 1) = {FF, FF, FF, FF, FF, FF, h}, (3 - 3 + 1) = {FF, FF, FF, -FF, -FF, -FF, h}, (3 + 2 + 2) = {FF, FF, FF, h, h, z, z}, 
(4 + 2 + 1) = {FF, FF, FF, FF, h, h, zj, (4 + 3) = {FF, FF, FF, FF, h, h, h} and (5 + 2) = {FF, FF, FF, FF, FF, h, h} cases. 
Again, first two are absolutely identical to similar cases in lower dimensions - isotropic case is stable as 
long as FF > 0 while (6 + 1) = {FF, FF, FF, FF, FF, FF, h} is unstable in F-direction. Of the remaining cases, 
(3 — 3 + 1) = {FF, FF, FF, —FF, —FF, —FF, h} is identical to (2 — 2 + 1) = {FF, FF, —FF, —FF, h} from (5+1)- 
dimensional case and is unstable for the same reasons, while (4 + 2 + 1) = {FF, FF, FF, FF, h, h, z} case is 
unstable in z direction since z E R from the solution (see [31]). Remaining three cases have nontrivial 
stability regions and now let us consider them. 

First case to consider, (3 + 2 + 2) = {FF, FF, FF, fi, F, z, zj, for the solutions of perturbed equations has 
6Hi = H5Hi making it stable iff FF <0. 
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Second of them, (5 + 2) = {H, H, H, H, H, h, h} case, has the following stability regions for perturbed 
equations: 


18aH‘^ - 1 
8aH 


> 0 


(9) 


and qualitatively they correspond to Fig. 1(a). Similar to previous cases with {{D — 2) + 2) splitting, this one 
also admits vacuum solution. Following the procedure from previous cases, one finds 4^^+16^^+40^+45 = 0, 
its solution ^ Ri —1.870 and after substitution to = —l/(8a(2^ + 3)) get a Ri 1/(5.9277^). Comparing it 
with (9), one can tell that vacuum solution is stable iff 77 > 0. 

Final case, (4 + 3) = {77,77,77,77, h, h, h} one, after solving perturbed equations demonstrates the fol¬ 
lowing inequality for the stability regions: 


1 10a77 ± 3 V24a2i72 _ q, 

2 a 


< 0 


( 10 ) 


where “±” corresponds to two branches according to the solution itself (see [31]). The shape of stability 
area resembles that from (3 + 3) = {77, 77, 77,/i, h, h} (6+)-dimensional case (see Eq. (8)) and so up to a 
factor represented in Fig. 1(b). The factor in question is - now the left and right wings are: for 77 < 0 we 
have a > 1/(2477^) while for 77 > 0 it is a > 9/(11677^). 

Again, for this splitting there also exists vacuum solution. Following the usual procedure, we have 
C^+6^^+11^^+12^+5 = Oand its solutions Ri —3.874 and ^2 ~ —0.743; having 77^ = —1/(477^(^^+6^+3)) 
we can derive ai ~ 1/(20.94477^) for the first and a 2 ~ 1/(3.62477^) for the second branches. This means, 
with above-mentioned bounds from (10), that the ^2 branch is stable for 77 > 0 while branch is stable 
only iff 77 < 0, which is quite similar to the (3 + 3) = {77, 77, 77, h, h, h} (6+)-dimensional case. 

This hnalize our study of the (7+l)-dimensional case. As expected, isotropic solution behaves “as 
usual” while the solutions with unconstrained exponents - (6 + 1) = {77,77,77,77,77,77, h}, (3 — 3 + 1) = 
{77,77,77, —77, —77, —77, h} and (4 + 2 + 1) = {77,77,77,77, h, h, z} - also behave as expected - they have un¬ 
defined perturbations in the unconstrained directions which makes them unstable. Of the remaining cases, 
(3 + 2 + 2) = {77,77,77, h, h, z, z} is stable as long as 77 < 0, (5 + 2) = {77,77,77,77,77, h, h} has “typical” 
behavior (see Fig. 1(a)) so as (4 + 3) = {77,77,77,77, h, h, h} - see Fig. 1(b). Vacuum solutions also behave 
similar to the (6+l)-dimensional case - vacuum (5 + 2) = {77, 77, 77, 77, 77,/i,/i} case is stable as long as 
77 > 0 while (4 + 3) = {77, 77, 77, 77, h, h, h} has two branches, one of them is stable while 77 > 0 and the 
second if stable iff 77 <0. 
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E. Conclusions 


This concludes our study of the Gauss-Bonnet solutions stability. In all cases we have two solutions 
which repeat themselves - isotropic solution and {{D — 1) + 1) splitting case. The former is stable iS H > 0 
while the latter is always unstable: “detached” dimension is unconstrained from the solution and through 
this the perturbations in that direction are unstable. 

In each dimension we have the case with stability regions looks alike each other (see see Fig. 1(a)): 
(2 + 2) = {H,H,h,h} case from (4-|-l) dimensions (see (4)), (3 + 2) = {H,H,H,h,h} case from (5-t-l) (see 
(6)), (4 + 2) = {H, H, H, H,h,h} case from (6+1) (see (7)) and (5 + 2) = {H, H, H, H, H,h,h} case from 
(7+1) (see (9)). One cannot miss that all of them have similar splitting - two dimensions are detached - 
((74 — 2) + 2) splitting. We discuss their properties more in the Discussions section. 

The same is true for (3 + 3) = {H, H, H, h, h, h} case from (6+1) dimensions (see Eq. (8)) and (4 + 3) = 
{H,H,H,H,h,h,h} case from (7+1) dimensions (see Eq. (10)). Again, both cases could be rewritten as 
((74 — 3) + 3) and look like presented in Eig. 1(b); we discuss them in detail in Discussions section. 

One more stable solution is (3 + 2 + 2) = {H, H, H,h,h, z, z} from (7+1) dimensions - it is stable iff 
77 < 0. 

Stability of vacuum solutions also have a pattern to follow through different dimensions - ((D — 2) + 2) 
vacuum solutions have only one branch thich is stable iff 77 > 0 while ((74 — 3) + 3) case has two branches 
- one of them is stable for 77 > 0 while the other one is stable as long as 77 < 0. Condition for stability of 
vacuum isotropic solutions is the same as for A-term solutions - it is 77 >0. 

All other cases are unstable for the same reason as ((D — 1) + 1) splitting case ~ in the solution we have 
one (or more) exponents unconstrained which leads to instability in that direction. 

III. SOLUTIONS WITH CUBIC LOVELOCK TERM 

This case is quite similar to the previous one. The full system of equations (1) is replaced with system 
which allows cubic Lovelock contribution: jth dynamical and constraint equations 


2 I + ^i) + S HiHk j + 8a I + Hf) H^Hi + 3 X) j + 


+144/3 (E(i^i +772) ^ 


_ _ IkHiHmHn + 5 Y. HkHiHmHnHpHr - A = 0; 

j} {k>l>m>n>p>r}i:j ) 


2 E 77,77, + 24a E + 720/3 E H.HjHkHiHmHn = A, 

i>j k>l>m>n i>j>k>l>m>n 


( 11 ) 
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where notations are the same as in (1) plus an additional cubic Lovelock coupling (5 is introduced and the 
procedure is exactly the same: we perturb the solution as Hi ^ Hi + 6Hi and Hi —)■ 6Hi and use the 
same stability criterium (2). Now let us consider (6+1)- and (7+l)-dimensional cases as they are two lowest 
possible dimensions with cubic Lovelock term. 


A. 6+1 

This case has solutions with eight different spatial splittings. It appears that first two cases - isotropic 
and (5 + 1) = {H, H, H, H, H, h} - are exactly the same as in Gauss-Bonnet case - exactly, isotropic one 
is stable as long as iL > 0 while the case with one detached dimension is unstable in that direction for the 
same reason as in the Gauss-Bonnet case. 

Apart from the last mentioned case, there are three more cases with unconstrained exponents - (4 + 1 + 
1) = {H, H, H, H, h, z}, (3 + 2 + 1) = {H, H, H, h, h, z} and (2 - 2 + 1 + 1) = {H, H, -H, -H, h, z}; all of 
them are unstable in these directions and (4 + 1 + 1) = {H, H, H, H, h, z} case is unstable in all directions. 

Remaining three cases have nontrivial stability regions in the {H, a, (3) space. Now we have three 
parameters and so representation of the stability regions is somewhat tricky. So we decided to plot areas in 
(a, (5) plane for some fixed iL - in that case the curves are subject to H. Generally it do not change the 
shape of the regions, only rescale and rotate them - we describe the behavior in each case. 

First case to consider is (4 + 2) = {H, H, H, H, h, h}, and the solution of the perturbed equations admits 

5Hi _ 1 UAI3H^ + UaH^ " ^ ^ n fl2) 

5Hi ~ 6 H{<opm +a) ^ ^ ^ 

We plot the corresponding stability regions in Fig. 2(a) and (b). In (a) panel we present the stability 
regions for H > 0 and H < 0 for a fixed H^ while in (b) we present the “evolution” of these regions with 
changing H^. From (12) one can see that the shape of the stability regions depend only on H^ while single 
H in denominator determines the sign - would it stable for iL > 0 or for H < 0. So in Fig. (c) shaded 
region corresponds to LI > 0 stability while white - to LI < 0. 

The second case with nontrivial stability area is (2 + 2 + 2) = {LI, H, h, h, z, z}; the stability condition 
for it yields 


6H, 9f3H{iaH^-l) 

6Hi a(18/ILl2 -a) ^ ’ 

and we plot them in Fig. 2(c) and (d). Again, in (c) panel we presented the structure if the stability regions 
while in (d) -- their variation with varying H. Similar to the previous case, shaded regions in (c) panel 
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Figure 2. Stability regions on the (a, /3) plane for solutions with cubic Lovelock terms. In (a) panel we presented 
typical behavior for (6+l)-dimensional (2 + 2 + 2) = {H, H, h, h, z, z} case (see Eq. (13)). Shaded region correspond 
to stability of H > 0 and white - to < 0. In (b) panel we presented the effect of varying H on the separation 
curves. In (c) panel we have typical behavior for {{D — 2) + 2) spatial splitting. This figure exactly corresponds 
to (4 + 2) = {H,H,H,H,h,h} case from (6+1) dimensions (see Eq. (12)), but up to a factor the situation with 
(5 + 2) = {H,H,H,H,H,h,h} case from (7+1) dimensions (see Eq. (15)) could also be described by this figure. 
Again, shaded region correspond to stability of iL > 0 and white - to H <0. Similar to the previous case, in (d) 
panel we presented the effect of H variation on the separation curves (see text for details). 

correspond to H > 0 stability and white - to H < 0. And the last case is (3 + 3) = {H, H, H, h, h, h} and 
this area is defined from 


6Hi 3 36PH^ - 2aH ± ^-72a^H^ + 1202^2 - - a 

iHi ~ ~2 + a 


(14) 


where “±” corresponds to two branches of the solution (see [31]). We presented the areas of stability in 
Fig. 3. In (a) and (b) panels we presented the structure of the stability regions - fine structure in (a) 
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Figure 3. Stability regions on the (a, /3) plane for (6+l)-dimensional solution with cubic Lovelock term and (3 + 3) = 
{H, H, H, h, h, h} spatial splitting. In (a) and (b) panels we presented the structure of the stability regions - fine 
structure in (a) and large-scale - in (b). Different shading correspond to stability to either H > 0 oi H <0; white 
region correspond to instability due to negativity of the radicand in (14). In (c) and (d) panels we presented the 
variation of separation curves with varying - again, in (c) we presented large-scale structure and in (d) - fine 
structure. 


and large-scale - in (b). Different shading correspond to stability to either H > 0 or H < 0; white region 
correspond to instability due to negativity of the radicand in (14). In (c) and (d) panels we presented the 
variation of separation curves with varying - again, in (c) we presented large-scale structure and in (d) 
- fine structure. 

This concludes our study of (6+l)-dimensional solutions stability with cubic Lovelock term. We described 
stability regions of three nontrivial cases - (4 + 2) = {H, H, H, H, h, h}, (2 + 2 + 2) = {H, H, h, h, z, z} and 
(3 + 3) = {H, H, H,h,h,h}. Isotropic case is stable as long as H > 0 - exactly similar to all previous 
cases. Finally, we detected instabilities due to unconstrained exponents for remaining four cases - (5 + 1) = 
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{H, H, H, H, H, h}, (4 + 1 + 1) = {H, H, H, H, h, z}, (3 + 2 + 1) = {H, H, H, h, h, z} and (2 - 2 + 1 + 1) = 
{H, H,-H,-H,h,z}. 

We also described the effect of varying on the geometry of the stability regions (as we plot stability 
regions on (a, j5) plane). One cannot miss substantial decrease in H > 0 stability with growth of (and 
overall stability in (3 + 3) = {H, H, H, h, h, h} case) - and this is true for all three cases considered. We 
address this point later in Discussions section. 


B. 7+1 

In (7+l)-dimensional cases we follow the procedure for stability regions representation from (6+1) di¬ 
mensions. This case is abundant with ten solutions (see [31]), but only three of them are stable with 
nontrivial stability regions. As usual, isotropic solution is stable as long as H >0. Six solutions - 
(6 + 1) = {H, H, H, H, H, H, h}, (5 + 1 + 1) = {H, H, H, H, H, h, z}, (4 + 2 + 1) = {H, H, H, H, h, h, z}, 
(3 + 3 + 1) = {ff, H, H, h, h, h, z}, (3 — 3 + 1) = {H, H, H, —H, —H, —H, h} and (2 — 2 + 1 + 1 +1) = 
{H, H, —H, —H, h, y, z} are unstable as one or more exponents are undefined which means that the pertur¬ 
bation in the corresponding direction is unconstrained which makes it unstable in that direction. Now let 
us consider the remaining three nontrivial cases. 

First of them is (5 + 2) = {H, H, H, H, H, h, h}; after solving perturbed equations we have the following 
constraint on the stability regions: 


5Hi _ 1 M8/3H^ + 16aH^ “ ^ ^ n flT 

One cannot miss familiarity between ( 12 ) and (15); if we plot regions from (15) we get the picture quite 
similar to Fig. 2(c). The same is true for the effect of H variation (see Fig. 2 (d)). 

Next case to consider is (3 + 2 + 2) = {H, H, H, h, h, z, z} and the resulting stability regions defined as 
follows: 


6Hi _ H{972(3^H^ - W8al3H^ - + 18/3) ^ . . 

Wi~ 324/32774 _ 36a^Fr2 + ^ ) 

We plot stability regions which correspond to the expression above in Figs. 4(a) and (b). One can see 
that it differs from (2 + 2 + 2) = {H, H,h,h, z, z} case in (6+1) dimensions (Eq. (13) and Fig. 2(a)). Shaded 
regions in Fig. 4(a) correspond to 77 > 0 stability and white - to 77 < 0. In (b) panel we depicted variation 
of separation curves with varying 77^. 

And the last case is (4 + 3) = {77, 77, 77, 77, h, h, h}; the solution of the perturbed equations holds 
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Figure 4. Stability regions on the (a, /3) plane for (7+l)-dimensional solutions with cubic Lovelock term. In (a) and 
(b) panels we presented the situation with (3 + 2 + 2) = {H, iJ, H, h, h, z, z} spatial splitting: in (a) panel we draw 
stability regions structure and in (b) - its variation with varying H. Shaded area in (a) panel correspond to H > 0 
stability and white - to iL < 0. Dashed line corresponds to discontinuity due to zero of the denominator in (16). 
In (c) and (d) panels we presented the situation with (4 + 3) = {H, H, H, H, h, h, h} spatial splitting: in (c) panel - 
stability regions and in (d) - their variation with varying H. Now in (c) panel shaded region corresponds to H > 0 
stability and double-shaded - to H <0. White region corresponds to instability due to negativity of the radicand in 
(17) (see text for details). 


6Hi _ 1 324(3H^ - lOaH ± 3 - 2l6al3H^ + 24a^H‘^ - - a 

Wi ~ ~2 + a 


(17) 


We plot the corresponding stability regions in Figs. 4(c) and (d). There in (c) panel we presented the 
structure of the stability regions. Shaded region correspond to H > t) stability while double-shaded - to 
H < t). White region correspond to instability due to negativity of the radicand in (17). In (d) panel we 
present the effect of H variation on the stability regions. 
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One can see familiarity and difference between this case and (3 + 3) = {H, H, H, h, h, h} (6+1)- 
dimensional: they both have instability region caused by negativity of the radicand in the solution, in 
both cases this region is defined by second-order curve, but the kind of curve is different in these cases. 
Also, stability for H < 0 case exist inside area between two crossing lines while for H > 0 solutions are 
stable on (almost in the (3 + 3) = {H, H, H, h, h, h} ( 6 +l)-dimensional case) the entire domain of definition. 

This finalize our study of the (7+l)-dimensional solutions stability. We report six out of ten solutions 
- (6 + 1) = {H, H, H, H, H, H, h}, (5 + 1 + 1) = {H, H, H, H, H, h, z}, (4 + 2 + 1) = {H, H, H, H, h, h, z}, 
(3 + 3 + 1) = {H, H, H, h, h, h, z}, (3 — 3 + 1) = {H, H, H, —H, —H, —H, h} and (2 — 2 + 1 + 1 +1) = 
{H, H, —H, —H, h,y, z} - to be unstable due to lack of constraint on certain exponents in correspond¬ 
ing solutions. Remaining four solutions are: isotropic solution is stable as long as H >0, (5 + 2) = 
{H,H,H,H,H,h,h} solution is stable if (15) is satished (corresponding structure of stability regions is 
similar to that in Figs. 2(c) and (d)), (3 + 2 + 2) = {H, H, H,h,h, z, z} case is stable while (16) is valid 
(see also Figs. 4(a) and (b)) and finally (4 + 3) = {H, H, H, H, h, h, h} case requires (17) for stability (with 
corresponding regions depicted in Figs. 4(c) and (d)). 

C. Conclusions 

This concludes our stability study of solutions with cubic Lovelock term taken into account. Unlike 
Gauss-Bonnet case now we have only two cases with known solutions (versus four for Gauss-Bonnet) which 
makes it harder to find familiarities. Still, some of them we have detected. 

First of all, two cases - isotropic and case with ((H — 1) + 1) splitting - behave exactly the same as 
in case of Gauss-Bonnet term - the former of them is stable iff 77 > 0 while latter is unstable due to 
unconstrained exponent in the detached dimension. The same reason - unconstrained exponent - leave the 
absolute majority of the solution found (see [31]) unstable. 

Three classes of spatial splitting - {{D — 2) + 2), ((Zl — 3) + 3) and ((77 — 4) + 2 + 2) - appear to be 
stable in both (6+1) and (7+1) dimensions. First of them depicted in Fig. 2(a) - actual figure is for (6+1)- 
dimensional case but (7+l)-dimensional fits it up to a factor. Second of them, ((71 — 3) + 3), drawn in Fig. 
3 for ( 6 +l)-dimensional case and Fig. 4(c) for (7+l)-dimensional. In this case we have different shape of 
separation curves, but some features still are similar, like, the forbidden region for both of them defined by 
second-order curve. Finally, ((71 — 4) + 2 + 2) case is also similar yet different - see Fig. 2(c) and Fig. 4(a). 

Two of these three cases have similar behavior of stability areas with changing 77: ((71 — 2) + 2) and 
((71 — 4) + 2 + 2) cases in both (6+1) and (7+1) dimensions have growing area of 77 > 0 and so shrinking 
area of 77 < 0 stability with increasing of 77^ (areas of 77 > 0 and 77 < 0 stabilities cover entire (a, P) 
plane and do not overlap). In contrast, ((71 — 3) + 3) case behaves differently - areas of 77 > 0 and 77 < 0 
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stabilities do not cover entire (a, fi) plane and do overlap. Additional features also differ in (6+1) and (7+1) 
dimensions ~ (3 + 3) = {77, H, H, h, h, h} (6+l)-dimensional case, presented in Fig. 3, have shrinking area 
of both H > 0 and H < 0 stabilities, while in (4 + 3) = {H, H, H, H, h, h, h} (7+l)-dimensional case area of 
H > 0 stability is expanding with growth of H^. 

Finally, vacuum counterpart of the solutions was not generally considered. Stability condition for vacuum 
isotropic solutions are the same as for A-term ones - they are stable as long as H >0. But from [31] one 
can clearly see that (4 + 2) = {H, H, H, H, h, h} and (3 + 3) = {H, H, H, h, h, h} in (6+1) dimensions as well 
as (4 + 3) = {H, H, H, H, h, h, h}, (3 + 2 + 2) = {H, H, H, h, h, z, z} and (5 + 2) = {H, H, H, H, H, h, h} in 
(7+1) dimensions also admit vacuum solutions. Their stability were not addressed due to high nonlinearity 
of the resulting equations. Indeed, the equations are obtainned as in the Gauss-Bonnet case, from A = 0 
condition from the constraint equation. In the Gauss-Bonnet case it results in bi-qubic equation, which 
could be dealt with. But in the cubic Lovelock case the constraint is reduced to bi-sextic equation which 
cannot be solved analytically any longer. It still could be solved numerically but in that case we would not 
be able to draw any generalizations so we decided to put this case aside for a while - we are going to return 
to it in one of the following papers. 


IV. DISCUSSIONS 

As we noted in the Introduction section, exponential solutions play important role among exact solutions 
in cosmology. In [29, 31] we reported exact exponential solutions in Gauss-Bonnet and Lovelock gravity and 
demonstrated that we found all possible solutions with varying volume. There appeared to be a number of 
different solutions and some of them could give rise to a successful compactification so additional investigation 
required to find out which of them are viable from different physical points of view. One of such tests to 
solutions is the stability check. So we perturb full Gauss-Bonnet (or more general Lovelock) system around 
solution to see how perturbations behave. With several parameters - {H, a) for Gauss-Bonnet and (77, a, /3) 
for solutions with cubic Lovelock term - for solutions which admit stability we can plot corresponding areas 
of stability on the parameters space. For Gauss-Bonnet case with only two parameters we simply use 2D 
(77, a) plane but in cubic Lovelock case there are three parameters and use of 3D plot could obscure the 
structure of stability regions. So we fixed 77 and plot regions in (a, /7) plane and put an additional plot 
which demonstrates the effect of changing 77 on the separation curves. 

Both Gauss-Bonnet and the case with cubic Lovelock term have two same features - first, isotropic case 
is always stable iff 77 > 0. This is true in any number if dimensions and with any order of highest Lovelock 
correction taken into account - corresponding solution of the perturbed equation is 5Hi{t) = with 

D being number of spatial dimensions. 
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The second feature is instability due to unconstrained exponent in the exponential solution. Indeed, 
in [29, 31] we reported a number of solutions with unconstrained one or more exponents, e.g., y E R. When 
one perturb such solution, the corresponding perturbation 5Hi{t) is excluded from the full perturbed system, 
leaving it unconstrained. We treat it as instability, as for stable solution we need damping perturbation and 
we cannot say this about unconstrained one. Again, formally we cannot call them “unstable” for the reason 
above, they would be rather in “neutral stability” class, but since we cannot call them “stable” either, we 
exclude them from stable solutions. So below we describe cases with nontrivial stability regions. 

It appears that in Gauss-Bonnet case there are two such cases - ((D — 2) + 2) and {{D — 3) + 3) spatial 
splittings. Typical stability regions for the former of them is presented in Fig. 1(a) and for the latter - in 
Fig. 1(b). By “typical” here we mean that in each particular dimension separation curves resemble this 
typical one up to a factor. Indeed, say, for {{D — 2) + 2) splitting Fig. 1(a) actually correspond to (4+1)- 
dimensional (2 + 2) = {H, H, h, h} case (see Eq. (4)), but one can check that other cases with {{D — 2) + 2) 
spatial splitting - (3 + 2) = {H, H, H, h, h} case from (5+1) (see (6)), (4 + 2) = {H, H, H, H, h, h} case from 
(6+1) (see (7)) and (5 + 2) = {H, H, H, H, H, h, h} case from (7+1) (see (9)) have expressions that meet (4) 
upto two factor which makes the corresponding figures resemble Fig. 1(a). 

For the (7+l)-dimensional Gauss-Bonnet case there is additional stable solution (3+2+2) = {H, H, H, h, h, z, z] 
- it is stable iff < 0 . 

In contrast to Gauss-Bonnet, the case with cubic Lovelock term has three spatial splitting which admit 
stability - ((T> — 4) + 2 + 2), ((L> — 2) + 2) and {{D — 3) + 3). But only one of them, namely, {{D — 2) + 2) 
has feature we explained in the Gauss-Bonnet case - resemblance of stability regions in different number of 
dimensions. Indeed, in Fig. 2(c) we presented (4 + 2) = {H, H, H, H,h,h} case from (6+1) dimensions (see 
Eq. (12)), but up to a factor the situation with (5 + 2) = {H, H, H, H, H, h, h} case from (7+1) dimensions 
(see Eq. (15)) could also be described by this figure. 

Two remaining cases - {{D — 4) + 2 + 2) and {{D — 3) + 3) ~ have different behavior in (6+1)- and 
(7+l)-dimensional cases. The former of them presented in Eig. 2(a) (( 6 +l)-dimensional (2 + 2 + 2) = 

{H, H,h,h, z, z} case (see Eq. (13))) and Eig. 4(a) ((7+l)-dimensional (3 + 2 + 2) = {H, H, H, h, h, z, z} 
case (see Eq. (16))). The latter is presented in Eig. 3 (( 6 +l)-dimensional (3 + 3) = {H, H, H, h, h, h} case 
(see Eq. (14))) and Eig. 4(c) ((7+l)-dimensional (4 + 3) = {H, H, H, H,h,h,h} case (see Eq. (17))). One 
can see that both of them have substantial differences between (6+1)- and (7+l)-dimensional cases. 

The last point to attend to is the variation of the stability regions with varying H in cubic Lovelock case. 
There we have familiarity between ((D — 4) + 2 + 2) and {{D — 2) + 2) cases - both of them have the entire 
(a, (3) plane covered but separated into two regions - with H > 0 stability and H < 0 stability regions. 
These regions, as said, cover the entire (a, /3) plane and do not overlap. With increasing the area of 
H > 0 stability gradually increase while H < 0 regions shrinks as presented in Eigs. 2(b), (d) and 4(d). 
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In contrast to two previous cases the —3)+3) case has quite different behavior. First, it has not entire 

(a, fi) plane covered - indeed, in both (6+1)- and (7+l)-dimensional cases we have radicand in the stability 
condition expression (see (14) and (17) correspondingly) whose negativity cause instability presented as 
white region in Figs. 3(a) and (b) and 4(c). Secondly, areas of 77 < 0 and H > 0 stability do overlap in that 
case - as areas of 77 < 0 and 77 > 0 stability presented as different shading in Figs. 3(a), (b) and 4(c) their 
overlapped region has double shading. And the final difference is the behavior of the separation curves with 
varying 77^: while (4 + 3) = {77, 77, 77, 77, h, h, h} case has “usual” growth of the 77 > 0 area with increase of 
77^ (see Fig. 4(d)), (3 + 3) = {77, 77, 77, h, h, K] case has shrinking 77 > 0 area with increase of 77^ (see Figs. 
3(c), (d)). This is the only case of all considered which has this property - all other cases have increasing 
77 > 0 area with growth of 77^. 

As we mentioned in the Introduction, of special interest are the solutions with three-dimensional isotropic 
subspace. In [29, 31] we reported a number of them, but current analysis demonstrate that not much of 
them are stable. In fact, in Einstein-Gauss-Bonnet case we can call only ((77 — 3) + 3) spatial splitting 
as stable (with 77 > 5); the same splitting is stable for the case with cubic Lovelock term. Formally, in 
Einstein-Gauss-Bonnet (7 + l)-dimensional case has solution with (3 + 2 + 2) = {77, 77, 77, h, h, z, z} spatial 
splitting and it is stable iff 77 < 0, but as we look for viable compactification scheme we want expanding 
three-dimensional subspace, not contracting, which makes this case inviable. Unlike Einstein-Gauss-Bonnet 
case, in the case with cubic Lovelock term (3 + 2 + 2) = {77,77,77, h, h, z, z} spatial splitting is quite stable, 
as a part of ((77 — 4) + 2 + 2) general case. So that among a number of reported solutions which could 
give rise to compactification only one solution in Einstein-Gauss-Bonnet case - ((77 — 3) + 3) - could be 
said to be stable. The same solution is stable in the cubic Lovelock case as well; additionally we have 
(3 + 2 + 2) = {77, 77, 77, h, h, z, z} solution in (7 + 1) dimensions. 

Vacuum solutions are obtained from A-term ones with A = 0 condition in the constraint which diminish 
the number of parameters by one. In the Einstein-Gauss-Bonnet case with A = 0 condition the constraint 
is reduced to bi-cubic equation and so this case is analyzed while in the cubic Lovelock theory it is reduced 
to bi-sextic equation and generally we cannot solve it analytically. So we presented analysis of the vacuum 
solutions stability in the Einstein-Gauss-Bonnet case but omitted them in the cubic Lovelock - we are going 
to return to this problem in the near future. 

To conclude, presence of cubic Lovelock term severely changes not only the abundance of the solutions 
(see [29, 31]), but also the situation with their stability. Indeed, the same spatial splittings in (6+1) and 
(7+1) dimensions have quite different stability conditions in Gauss-Bonnet and case with cubic Lovelock 
term taken into account. We can also conclude that we described truthfully Gauss-Bonnet case and have all 
grounds to generalize solutions found on any number of dimensions, while for cubic Lovelock term case there 
are still differences between different number of dimensions; sometimes these differences are quite severe, as 
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in ((li) — 3) + 3) splitting case. So perhaps we need additional study of higher-dimensional cosmologies in 
cubic Lovelock case to draw any generalized conclusions. 
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